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I. 



INTRODUCTION 



Inventory control is a pivotal activity of any logistics 
organization. Multi-item inventory systems encompass trade- 
offs in balancing customer service needs with operating 
costs. This management task is particularly challenging 
in the military setting where item availability often 
affects mission readiness. Clearly, one objective of any 
inventory doctrine is to succeed in making those decisions 
which minimize operating costs while providing an acceptable 
level of service for a forecasted rate of demand. The 
demand forecasting process and its influence on inventory 
system effectiveness are the subjects of this study. 

First, some popular forecasting methods are profiled in - 
terms of six key evaluation measures. 

In recent decades, a wide variety of forecasting methods 
has emerged. Generally, they can be assigned to one of two 
taxonomies: qualitative or quantitative. Qualitative 

techniques are regarded as the more subjective of forecasting 
approaches. Usually conducted in a setting where historical 
data is unavailable, this class often employs expert opinion 
in constructing a forecast. The Navy Supply System uses 
such a method, called Best Replacement Factors, when esti- 
mating initial stock levels for a new item of inventory. 

In contrast, quantitative methods make extensive use of 
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historical data. There the data serves as input for various 
types of mathematical models which compute the required 
forecast. Not surprisingly, advances in computer technology 
tend to popularize the quanti tative-oriented forecasting 
methods. Rather than spending considerable time synopsizing 
the more common forecasting techniques, the interested 
reader is invited to consult the existing literature for 
supporting detail (see Makridakis and Wheelwright [Ref. 1]). 
However, it is interesting to classify a few of the widely 
used techniques in terms of their cost, accuracy, type, 
applicability, data pattern and time horizon characteristics. 
Adapted from Wheelwright and Makridakis [Ref. 2] , the pro- 
file provided by Tables 1 and 2 assists in selecting the 
most appropriate method for a given forecasting requirement. 
For example, the widespread appeal, among industry and the 
military, for the exponential smoothing technique is 
apparent. To utilize exponential smoothing a manager need 
have only three data elements: the most recent observation, 

the most recent forecast and a weighting parameter. This 
data storage consideration has been of primary importance 
to multi-item inventory systems where demand forecasts are 
routinely prepared for several thousand items. Combining 
such features as low data processing and storage costs to- 
gether with high applicability, exponential smoothing appears 
as a rational choice for the Navy’s forecasting method. In 
contrast, the Box-Jenkins technique is not suitable for 
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UICP use. Although it offers great accuracy, the Box- 
Jenkins method's excessive data processing and data storage 
requirements render its costs prohibitive. 

However profiled in terms of these six factors, the 
best measure of effectiveness for a forecasting method is 
the economic benefit it provides to its dependent inventory 
control system. The identification of a method which incurs 
low data processing costs while contributing to minimal holding 
and stockout costs remains the goal of the inventory fore- 
caster. This notion will be pursued further. 

Before examining the cost effects of forecast accuracy 
several issues will be analyzed. The following chapter re- 
views the forecasting model currently installed at one of 
two Navy Inventory Control Points, the Navy Ships Parts 
Control Center (SPCC) , and also introduces an alternate 
computer-oriented model called, "Focus Forecasting." Chapter 
III reveals that items having erratic demand patterns con- 
stitute an appreciable portion of the ICP's inventory 
population. Chapters IV and V identify both the cost 
effects and some sources of extreme variability in the UICP 
demands. Lastly, Chapter VI uses a computer simulation to 
examine the cost effects of forecast accuracy. 
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II. TWO FORECASTING MODELS 



A. EXPONENTIAL SMOOTHING 

As previously indicated, exponential smoothing satis- 
fies the requirements for a forecasting procedure which 
demands low computer storage and run time. Accordingly, 
the UICP model utilizes a refined version of exponential 
smoothing in establishing, among other things, the parameters 
for the leadtime demand distributions. Employing a modifi- 
cation of the continuous-review inventory formulation found 
in Hadley and Whitin [Ref. 3] , the UICP model applies the 
forecasted leadtime demand in determining both economic 
order quantities and reorder points. 

Exponential smoothing methods, originally advocated by 
Brown [Ref. 4], are a geometrically weighted sum of all past 
demands with the greatest weight applied to the most recent 
observation. Mathematically a forecast is calculated as: 



NEW FORECAST = LAST FORECAST + WEIGHT (LAST OBSERVATION 

- LAST FORECAST) 

Notationally this may be stated as 



F 



t+1 



F t + o(D fc - F t ) ; t > 1 



which reduces to 
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F 



t+1 



aD 



t 



+ (l-a)F 



t 



t > 1 



Here a is known as the smoothing weight and is normally 
assigned a value between zero and one. The accuracy of the 
exponentially smoothed forecast depends strongly on the 
chosen a value. In normal practice, the a smoothing con- 
stant is either selected arbitrarily or suggested through 
exhaustive sensitivity analyses. With D^_ again denoting 
the actual demand observation recorded in period t, the 
above equation can be represented in recursive form as: 



From this recursive form it is easily seen how exponential 
smoothing dilutes the effect of the older observations. 

The data-processing simplicity of this technique is evident 
since one stored value, F fc+ ^ , replaces the entire block of 
t demands. One point of concern can be the initial or seed 
value for F fc . Since a F value is needed before the next 
F fc+ ^ forecast is prepared, a seed value must be found to 
initialize the process when t = 0 . Several solutions to the 
seed selection problem are available. Some of the simpler 
solutions recommend using the first observation as a seed or 
dividing the time series data into two parts with the first 
part reserved for initial estimation purposes (average, 
least squares estimation, etc.). In practice the seed problem 



F 



t+1 
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is only of theoretical concern. (See Makridakis and 
Wheelwright [Ref. 1].) Generally, the forecasting process 
will be in operation long enough to suppress any dependency 
on the seed. 

The simple form of exponential smoothing described above 
is essentially an approximation of a moving average fore- 
casting process. An inherent drawback of this model is that 
it is relatively insensitive to recent trend changes. This 
weakness, left uncorrected, would frequently result in 
biased forecasts. Recognizing this limitation, the UICP 
model refines the process by incorporating two types of 
demand filters. (See Basic Inventory Manager's Manual [Ref. 
5].) First a trend test is conducted to detect sustained 
(>_ 3) changes from the underlying pattern. The trend 
statistic, TR, is a ratio test consisting of: 

2 (SUM OF LAST TWO OBSERVATIONS) 
iK (SUM OF LAST FOUR OBSERVATIONS) 

A trend is considered present when either: 

TR > 1.1 and (D fc >_ F fc , D t+1 1 F t+1* 

or 

TR < 0.9 and CD fc < F fc , D t+1 < F t+1 ) 
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When either of these trend conditions exist, the smoothing 
constant is modified from its usual .10 value to a new, 
"heavier" weight of .30. This causes the next forecast to 
be more directly influenced by the most recent observation. 
Next, a second filter is used to check for outlier observa- 
tions. This filter computes a control tolerance band around 
using multiples of mean absolute deviation (MAD) , the 
average of the absolute difference between actual and fore- 
casted demand. For further discussion on the exact UICP 
use of MAD see the Basic Inventory Manager's Manual [Ref. 

5] . If the most recent demand observation lies within a 
band of width 7.5 MAD about F fc , the process is considered in 
control and no modifications are necessary. If a single 
out-of-control condition has been indicated by the filter, 
the outlier demand is ignored and the forecast is left 
unchanged, i.e., F fc+ j_ = Further, if two successive 

demand observations lie on the same "side" (high or low) 
of the tolerance band, the next forecast is calculated as 
the average of the two cohort outlier demands, i.e., 

F fc+ ^ = . 5(D t _^ + D fc ) . This condition is known as a "step 
increase/decrease . " 

This completes the discussion of the SPCC exponential 
smoothing model. Next, "Focus Forecasting" is introduced. 

B . FOCUS FORECASTING 

Focus forecasting is a new forecasting approach first 
advocated by Smith [Ref. 6] . His concept requires the 
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dynamic simulation capabilities of modern computers in 
preparing each forecast. The methodology is somewhat 
inviting due to its overriding simplicity. Focus fore- 
casting emphasizes a straightforward, flexible design in 
acquiring user understanding and confidence. The require- 
ment for transparent forecasts which are derived from simple 
strategies and which capitalize on recent advances in com- 
puter technology motivated Smith in creating the concept. 

The mechanics of the process consist of four operations: 
backforecasting , selection, application and repetition. 

First, employing a dynamic evaluation routine, the computer 
identifies from a corpus of simple strategies the one which 
would have best forecasted the preceding period's demand. 

Next, this selected strategy is used in preparing the up- 
coming period's demand forecast. Lastly, the process is 
repeated until "optimal" strategies have been identified and 
applied for each inventory item. Recent strides in computer 
technology make possible Smith's procedure which, when imple- 
mented for large inventory systems , requires great internal 
processing speed. For clarity, a conceptual example of the 
focus forecasting algorithm is presented. 

• Item Demand History: 

Time Period (Quarters) :12345678 

Units Demanded (D ) : 596 388 527 259 270 363 357 250 

Required: F fc where t = 9. 
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Strategy Corpus: 



STRAT(l) — "LAST PERIOD THIS YEAR" 



F = D 
t t-1 



STRAT ( 2) — "AHEAD/BEHIND THIS QUARTER LAST YEAR" 



F = (D n x D . ) /D 
t t-1 t-4 ' ' t-5 



Routine: 



Step 1 — Backforecast for t-1, i.e., 
STRAT (1) = Fg = D ? 

= 357 



STRAT (2) 



Fg = (D 7 xD 4 )/D 3 

= (357 x 259 ) /527 

= 175 



compute F 



t-1 



Step 2 — Select item — "optimal" strategy 



DIFF(l) 



DIFF (2) 



ABS ( Dg - STRAT (1)) 
ABS (250 - 357) 

107 

ABS (Dg - STRAT ( 2 ) ) 
ABS ( 250-175) 

75 



since 75 < 107 item "optimal" strategy is STRAT(2). 



Step 3 — Apply selected strategy 



F t = (D t-l XD t-4 )/D t-5 
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F 9 " D 8 X V D 4 

= (250 x 270 ) /259 

= 261 



Step 4--Repeat for next item 
(go to step 1) 

Smith has implemented a seven-strategy version for a major 
commercial wholesaler of hardware. There the system pre- 
pares demand forecasts for over 100,000 items each month. 
Smith maintains that his procedure of adapting a series of 
forecasting approaches to item demand will significantly 
outperform a single-formula process such as the UICP exponen- 
tial smoothing model. Before this claim is examined, we 
briefly outline how focus forecasting might be successfully 
applied to military use — specifically leadtime demand 
parameterization . 

As an initial, albeit crude, attempt at implementation, 
six simple strategies derived from three general forecasting 
categories are to be used. The first two strategies come 
from the so-called, "Naive" class. These are the previously 
described "last period this year" method and its cohort, 

"this period last year." (Notationally — STRAT(l) = D 

STRAT(2) = D t _^.) The first method acknowledges trend; the 
second seasonality. The next three strategies are selected 
from the class of moving averages, all of which model hori- 
zontal demand patterns. The moving averages will be computed 
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as 2, 4 and 8 period averages. (Notationally — STRAT (3) = 

• 5(D t-l + °t-2 } ; STRAT ( 4 ) = • 25 ( D t _ 1 + ... + d 4 ) ; 

STRAT (5) = .125(D t _^ + ... + D t _g).) Finally, as proposed 
by Bates and Granger [Ref . 7] , a composite forecast is 
adapted in order to improve forecast accuracy by capturing 
information from each forecast strategy. (Notationally — 
STRAT (6) = .2 (STRAT (1) + ... + STRAT ( 5) ) ) . 

The next chapter compares the UICP exponential smoothing 
and focus forecasting models using empirical data. 
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III. MODEL COMPARISON 



A. EMPIRICAL DATA 

To evaluate the two forecasting models under considera- 
tion, a nine year demand history was retrieved. The Opera- 
tions Analysis Department at the Navy Fleet Material Support 
Office (FMSO) followed a stratified sampling procedure in 
preparing a random sample of 522 repairable and 4530 con- 
sumable items. Each of the 5052 items was represented by 
one master data record followed by several subrecords. The 
master record contains descriptive information such as 
national stock number, replacement price, etc. Each sub- 
record contains demand quantity and demand Julian date 
information. The complete record format can be found in 
[Ref. 8] . For convenience, sequential dates ranging from 
0001 to 3285 replaced the Julian dates for the nine year 
period. A FORTRAN computer program aggregated the subrecord 
demand data into 36 quarterly "buckets." (A complete listing 
of the principal FORTRAN IV source codes used in this thesis 
is available in Appendix A) . Negative demand quantities 
were assumed to be the result of customer cancellations. 

When negative demands were encountered, the quarterly demand 
balance was appropriately reduced. When adjusting for such 
cancellations, however, the quarterly balances were forced 
to maintain strictly non-negative values. 
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B. EVALUATION CRITERIA 



One aim of a forecasting process is the minimization of 
the total forecast error incurred over time. When stochas- 
tic demand rates are involved, a natural assumption is that 
superior inventory control requires very accurate forecasting. 
This assumption seems reasonable since, in the extreme case 
where demand rates are deterministic, existing inventory 
models achieve zero stockout cost and optimal ordering and 
holding costs. (See [Ref. 3].) Unfortunately, forecast 
accuracy lacks an absolute standard of measurement. In- 
stead, the decision-maker is free to choose from a wide 
variety of evaluation schemes. Conceivably, each scheme 
could rank the forecasting models differently. This section 
briefly describes the evaluation criterion selected for 
comparing the focus forecasting and SPCC exponential smoothing 
models. Chapter VI presents an alternate effectiveness 
measure — the cost impact of forecast error. 

Traditionally, squared forecast errors (SFE) have been 
a useful determinant of accuracy with mean squared error 
(MSE) serving as a popular choice. MSE is defined mathe- 
matically as: 



MSE = | l ( F t -DJ 2 

n t=l r r 

MSE is functionally related to variance and also enjoys 
wide acceptance as a measure of "closeness." Squaring a 
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forecast error provides two advantages. First, the alge- 
braic sign of the error is disregarded. This avoids distor- 
tion caused by offsetting positive and negative errors. 

Second, the squaring operation penalizes large forecast 
errors . 

SFE measures are not without drawbacks. For example, 
what is the preference ordering for two models posting 
MSE's of 81.49 and 75.19? The second forecasting model 
appears more accurate but with what significance? Another 
drawback of MSE is that it does not facilitate comparison 
across different time intervals. These limitations not- 
withstanding, three SFE measures will be used: mean, standard 

th 

deviation and the 9C) quantile. 

The next effectiveness criterion is designed to check 
for biased forecasts. Bias distorts forecast accuracy 
through systematic overforecasting or underforecasting. 

Several statistics drawn from the distribution of forecast 
errors help identify an unbiased forecasting model. Any 
reasonably unbiased model will exhibit two attributes. 

First for large samples, the forecast errors should have a 
mean of zero. Second, as an indication of symmetry the 
median forecast error should approximate its mean. The 
standard deviation of forecast errors will also be included 
as a measure of dispersion. 

A third important consideration in determining forecast 
accuracy is how well the model performs for the higher 
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priced inventory. The manager may benefit little from a 
model providing low bias and dispersion for only the lower 
priced items. Therefore, forecast errors will be weighted 
by item replacement price (PWFE) and tested by the three 
bias measures described above. 

Finally, the correlation between forecasted and actual 
demands is considered. Often two time series are closely 
related but not in a statistically dependent sense. The 
correlation coefficient can measure the strength of the 
linear relationship between two random variables. Also 
known as Pearson's product moment coefficient and denoted 
by r, the correlation coefficient is defined as: 

n 

l (D-D) (F -F) 

_ t=l _ 

< Z < D -D, 2 Z (F t -F) 2 )^ 2 
t=l t=l 

where D is mean actual demand; F is mean forecasted demand ; 
r e [-1,1] . 

A forecasting mechanism which parallels changes in the 
demand pattern is highly desirable. A high r value would 
typically indicate such behavior although it must be re- 
garded with caution. For example, a high r value may 
indicate a spurious relationship caused by chance or by 
the elimination of a third explanatory variable. Nonethe- 
less, r remains a recognized measure of the tracking capa- 
bility of a forecasting model. 
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In summary, four classes of evaluation criteria will 
be reported: 

1) SFE — (MSE, SDjP^q ) as measures of closeness, 

2) FE — (mean , median , SD) as measures of bias, 

3) PWFE — (mean , median, SD) as measures of price weighted 

bias , 

4) CORRE — (r) as a measure of association and direction. 

Before reporting the results of the model comparison, 

some amplifying remarks are offered. The SPCC exponential 
smoothing model was translated into the ES FORTRAN source 
code found in Appendix A. This code contains all of the 
provisions described in Chapter II. The seed values were 
computed as the average of the first four quarters of demand. 
To dilute the effect of this selection, the program was 
stabilized over demand years two and three. The demands 
representing years four through seven were examined in terms 
of the selected effectiveness measures. The FOCUS FORTRAN 
source code represents the previously described focus fore- 
casting model. For the six simple forecasting strategies 
selected, neither initialization nor stabilization was re- 
quired. To preserve an equal footing, however, demand years 
four through seven were again used for comparison to the 
SPCC model. Note that each strategy must require no more 
than eight quarters of demand history. This constraint is 
imposed since UICP data files presently access only two 
years of demand history. Also note that the FOCUS program 
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lacks a high demand filter. Some basic filters are offered 
by Smith but none were coded for the initial trial run. 

The tables which follow show the results of testing the 
SPCC exponential smoothing and focus forecasting models 
using four years of actual UICP demand history. To facili- 
tate model comparison between repairable and consumable 
inventory items, the demand history was separated into two 
populations. These tables record each forecasting model’s 
performance in terms of the previously described forecasting 
effectiveness measures. 

C. FINDINGS 

The results in Tables 3 and 4 do not suggest a clear 
modelling preference. However, they do reveal some unexpected 
results which led to a redirection of this study. First 
the MSE and SFE Pgg values indicate that the presence of 
large outliers severely distorts MSE. The use of Pg^ as 
a more stable measure does not indicate a significant dif- 
ference in model closeness. Second, both models succeed 
in providing unbiased forecasts with the focus forecasting 
model generally producing a smaller price weighted forecast 
error. Next, there does not appear to be any significant 
difference in the model's tracking ability as measured by 
r. However, both models appear more successful at tracking 
consumables than repairables. 
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Repairable Items — Sample Size 2088 (522 x 4 Quarters) 
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Correlation 



Consumable Items — Sample Size 18120 (4530 
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Correlation 



The next originally planned task was to fine tune focus 
forecasting using high demand filters and inject alternate 
forecasting strategies for subsequent trial runs. However, 
this task was interrupted when demand patterns , such as the 
three actual demand histories shown below, were observed. 



DEMAND HISTORY 

QUARTER 



ITEM 


1 


2 


3 


4 


5 


6 


7 


8 


1 


76 


5843 


18798 


15 


58 


746 


19 


0 


2 


6 


2 


6 


78 


0 


8 


10 


3 


3 


15 


44 


64 


7 


100 


604 


130 


2239 


It was 


startling to 


observe a 


significant 


amount 


of 


varia- 



bility in both the repairable and consumable demand histories. 
To gain a rough estimate of the magnitude of this dispersion 
the average coefficient of variation (STANDARD DEVIATION [D] / 
EXPECTED VALUE [D] ) was computed. The two data sets 
recorded average coefficients of: 



Repairable 


COEFF 


VAR = 


3.80 


Consumable 


COEFF 


VAR = 


3.58 



Since coefficients of variation greater than 1.0 are usually 
regarded as an indication of excessive dispersion [Ref. 1J , 
it appears that the demand histories profiled above are the 
rule rather than the exception. Recognition of the data's 
extreme variability is further suggested by SPCC's 1980 
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change to their model's high demand filter. The control 
tolerance band was revised from its previous F ± 3.75 MAD 
to a widened range of F fc ± 7.50 MAD. These erratic demand 
patterns challenge the foundation of those inventory models 
(such as SPCC's) which assume a stationary demand rate from 
quarter to quarter. Thus, two new research directions were 
suggested: (1) identify the causes of erratic UICP demand 

and (2) assess the cost impact of forecast inaccuracy on the 
type of inventory models implemented by the ICP's. 
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IV. THE ECONOMIC BENEFITS FROM VARIANCE REDUCTION 



Before studying the causes of erratic UICP demands, it 
is proper to first identify the marginal benefits which can 
be derived from reducing demand variance. 

The stochastic backorders model developed in [Ref. 3] 
contains three cost elements: order, holding and stockout 

costs. Mathematically these are represented in the total 
annual variable cost equation as: 

oo 

K(Q,r) = AA/Q + IC(Q/2 + r-y) + irA/Q[ / xh(x)dx - rH(r)], 

r 

or 



TVC = ORDER + HOLDING + STOCKOUT COSTS 

For our purposes, K(Q,r) is assumed to be both differen- 
tiable and jointly convex in Q and r. Under these assumptions, 
the values Q* and r* which minimize K(Q,r) are determined 
by the methods described below. It should be noted, however, 
that joint convexity depends on the particular distribution 
of h(x) . Brooks and Lu [Ref. 9] and Veinott [Ref. 10] 
show that K(Q,r) can be nonconvex when h(x) is non-decreasing 
and r < y. (Nonconvexity may lead to a failure in the 
optimization technique which follows.) 
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Under joint convexity the optimal values Q* , r* (where 
0 < Q* < «>, 0 < r* < °°) must satisfy the equations: 



3K A A , — , \ \ , IC 

_ = o = -j(A + nn(r>) + -j- 



}* (r) = [2X(A + un(rJ/IC] 1/2 



H ■ « 



= ^[|^-( / (x-r) h(x) dx) ] + IC 



= -^j-[-H(r)] + IC (using Leibniz's Rule) 



-*■ H (r* (Q) ) = QIC/ttX 



(note it must be sufficiently large such that H(r*(Q)) < 1) . 

Computing a numerical solution for Q*(r) and r*(Q) poses a 
problem of composite dependency. That is to solve for Q* 
we need to know n(r) and thus r. Secondly, to solve for r* 
requires a knowledge of Q. As a practical and accurate 
alternative, a five-step numerical iteration routine is 
invoked. 

1) Check for a unique solution, i.e., if 

[21 (A+Tty ) /IC] < ttA/IC unique Q*,r* 

2) Assume n(r) = 0 and starting with i = 1, 

compute = [2AA/IC]^ /2 
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3) 



Compute r i from H(r^(Q i )) = Q^IC/ttX by consulting 

"tail" distribution tables of H(x) 

4) Compute by using r^ to find n(r^) , i.e., 

compute = [2X (A+irn (r^) ) /IC] 

5) Repeat steps 3) and 4) until ,r ±-hl conver 9 e on Qj_,r^. 

2 

To determine the effect of demand variance (a ) on K 
requires a knowledge of 3K/3a where 

K(Q,r,a) = XA/Q(r,a) + IC(Q(r,a)/2 + r(a)-y) + TrXn(r,a)/Q(r,a) 

Unfortunately, developing a determinant form of 3K/3a is 
more complicated than the earlier task of solving for Q* and 
r*. The complication stems from the composite dependency 
among a, Q and r. That is, 3K/3c depends on 3n(r)/3o which 
depends on r which in turn depends on Q. For example, if 
the leadtime demand distribution is assumed normal with 
mean y and standard deviation a then 

3n (r) 3 r r°° / x . , . , , 

— ~ = do [ > (x-r) <J) (x , a) dx] 

r 

/ - §£ <J) (x, a) dx + / (x-r) do 

r r 

2 2 

where <J>(x,a) = ( l//2Tra) exp [ (x-y) /2a ] and r is the solution 
to H(r,a) = ICQ(r,a)/TrX. 

To remedy this complication and to graphically illustrate 
the relation between K and a, a sensitivity analysis was 
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conducted using the FORTRAN program entitled INVENTOR. For 
a set of input parameters which remained fixed throughout 
the analysis 1000 realizations of o were used in computing 
the three cost elements of K. The particular set of input 
parameters (A = 500, \ - 180, IC = 15, tt = 500, y = 90) 
was selected to ensure joint convexity in Q and r. The 
1000 values of a are equally spaced over the interval [10,350] . 
The upper and lower limits of this interval were selected to 
reflect the coefficients of variation found in the UICP 
sample. The upper limit of this range was specified to 
coincide with the average coefficient of variation (approxi- 
mately 3.80) reported in Chapter III (i.e., 350 ~ 3.80x90). 
Graphical representations of the sensitivity analysis are 
found in Figure 1-4. To facilitate comparison of the changes 
in each of the four cost categories, each axis shown in 
Figures 1-4 was normalized in the following manner: 

(1) Ordering, holding, stockout and total variable costs 
were each computed 1000 times for the 1000 different 
values of a. 

(2) The first value of each cost (where a = 10) was used 
as a basis value. 

(3) Each of the four sets of the 1000 cost values was 
then divided by its respective basis to produce the 
vertical coordinates used in the graphs. The 1000 a 
values were normalized by dividing each by a = 10 
thereby producing the horizontal coordinates found 
in the graphs . 
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NORMALIZED CHANGES IN ORDERING COSTS 




Figure 1. Ordering Cost Sensitivity to Sigma 
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Figure 2. Holding Cost Sensitivity to Sigma 
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Figure 3. Stockout Cost Sensitivity to Sigma 
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IN TOTAL 




Figure 4 . Total Cost Sensitivity to Sigma 
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These graphs highlight, in a general sense, the theoreti- 
cal cost effects of variability in leadtime demand. Figure 
1 shows that as a increases over the relevant range the 
ordering costs developed in the Hadley -Whi tin <Q,r> model 
decay in a nonlinear manner. In contrast. Figures 2 and 3 
show the corresponding effects on the holding cost and 
stockout cost components of the <Q,r> model. There both 
cost elements increase in a roughly linear fashion over the 
relevant range with (3BC/3a) > k x (3HC/3a) (where k = 4.5) . 

The reaction of total variable costs to a has been included 
to show the overall effect. It should be remembered that 
the effects of each of the other three cost elements will 
be dampened when translated into total variable cost. 

This sensitivity analysis indicates that a study of the 
sources of leadtime demand variance would help reduce holding 
costs but would also generate particularly beneficial effects 
on stockout costs. The identification of the causes of 
demand variance is the subject of the following chapter. 
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V. HYPOTHESIZED CAUSES OF DEMAND VARIABILITY 



Demand randomness results from two sources. There is 
the truly patternless variation that is unpredictable and 
uncontrollable, but there is also the variation which re- 
sults from the way data are collected or reported during a 
given operational schedule. This second source of variation 
is the subject of analysis of this chapter. 

Erratic items constitute an appreciable portion of the 
military's inventory population. Silver [Ref. 11] addressed 
this issue and it is also demonstrated by the coefficients 
of variation found in the UICP sample data discussed above. 
Silver maintains that imperfections in the wholesale manage- 
ment information reporting system are a primary cause of 
erratic demand. (FMSO calls their inventory reporting network 
the Transaction Item Reporting (TIR) System.) 

The TIR system is a tri-level organization which sub- 
scribes to a policy of geographical responsibility for supply 
support. The first level of the TIR's multi-echelon network 
is the end-user or consumer level (e.g., a shipboard supply 
department) . Level two is composed of several retail out- 
lets such as the at-sea Mobile Logistics Support Force 
(MLSF) ships and the shore-based stock points (Naval Supply 
Centers, industrial rework facilities. Naval Air Stations, 
etc.). The final tier contains the Navy's two wholesale 
activities — the Inventory Control Points. 
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Under the most common scenario, a shipboard supply 
department receives a routine demand request from one of 
its customer departments. If the demand reduces the ship's 
available stock past the reorder point or if the demand 
is for an unsupported item, a replenishment requisition is 
submitted to the nearest retail outlet by one of two means. 

For ships operating at sea with an assigned MLSF unit, 
the replenishment requisition is submitted directly to the 
MLSF ship. For ships operating without MLSF support or for 
those which are inport when the demand occurs, the replenish- 
ment requisition is deposited with the nearest shore-based 
retail outlet. Clearly for situations where a customer-ship 
operates at sea for extended periods without MLSF support, 
the replenishment requisition encounters substantial delays 
before reaching the TIR's retail level. Under either case 
when the requisition ultimately reaches a retail activity, 
the TIR system activates and provides transaction status to 
all levels. When the available stock at the retail outlet 
reaches its reorder point, a collective replenishment requi- 
sition is submitted to the ICP. This requisition reflects 
the current and anticipated future demands from the retail 
outlet's geographically assigned customers. In this cas- 
cading fashion every demand which originates at the consumer 
level ultimately reaches the wholesale level where the ICP 
is tasked with the system-wide reorder and budgetary 
responsibilities . 
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A simulation was conducted to illustrate how the TIR 



organizational structure may cause truly non-erratic con- 
sumer level demands to be observed as the highly erratic 
demand actually seen at the ICP level. (See Lewis and Uribe 
[Ref. 12] and FORTRAN source code, DEMSIM.) The framework 
for this simulation was adopted from Shields [Ref. 13] . 

The simulation modified the present TIR demand reporting 
scheme by inserting an artificial reporting channel which 
reports all consumer demands directly to the ICP as they 
occur. To date this reporting medium has been reserved for 
only high priority requisitions. This alternate reporting 
channel provides two valuable capabilities. First by directly 
reporting user-demands, the ICP's records will accurately 
reflect the demand patterns occurring on the consumer level. 
Second, the ICP recording operation is expedited by elimin- 
ating the reporting delays incurred when ships operate with- 
out MLSF services. Together these accuracy and timeliness 
considerations can have a big impact on the effectiveness 
of an ICP. 

The structure for the DEMSIM program is as follows: 

1) assume a 30 ship population with each ship operating 
independently over an 18 year active service life. 

2) initial ship inport and at-sea assignments were 
randomly selected. (The terms "inport" and "at-sea" 
are meant only to describe when a retail supply 
activity is readily available) . 
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3) assume inport periods follow a Uniform [5,30] dis- 
tribution and at-sea periods follow a Uniform [5,45] 
distribution. 

4} generate demands in a random manner following a 

compound Poisson process where demand interarrival 
times are Exponential and demand requisition sizes 
are Poisson. 

5) record demands according to: 

a) the indirect reporting system presently used 
by ICP. 

b) the artificial system where the ICP receives 
observations directly from the consumer. 

6} aggregate each demand record into quarterly groupings 
and compute the variance of each time series . 

7} repeat the process 50 times. 

The result of the simulation supports a similar finding 
to that of Shields [Ref. 13] , that demand variance is sub- 
stantially reduced using the alternate reporting method. 
After 50 replications the direct reporting scheme demon- 
strated an average variance of 3050 as opposed to 5560 for 
the indirect ICP reporting system (.a 45% reduction] . 

This simulation is meant to highlight only one cause of 
erratic demand. Presumably other changes in the demand 
recording system could yield additional improvements . 
Examples of such changes are: 

* Limited Demand History-continued strides in computer 
technology will make storage of longer demand histories 
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an economic reality. Greater depth in the demand 
history would facilitate aggregation by other than a 
quarterly basis. A change in the base time period 
could provide a smoother stream of observations. 

* Customer Identity — The interaction of large customers 
(e.g., naval bases) who order more product less often 
with small customers (e.g., destroyers) who order 
less product more often might cause high variability. 

It may prove advisable to monitor a few large 
customers on an individual basis. 

* Customer Buying Habits — A practice of batch ordering 
could easily be a contributing factor. Actual consump- 
tion of material may follow a stable pattern, however, 
due to such factors as funding and increased operational 
commitments material is often ordered in large lots. 

* Product Identity — The disaggregation of products 
into their life cycle states (growth, steady state, 
gradual phase-out) might produce different and dis- 
tinct demand patterns . 

Whether such changes could or even should be implemented 
largely depends on the results gleaned from comparing the 
inventory cost savings against the cost of altering the 
demand recording system. This study does not undertake 
such a task but focuses instead on a more pertinent issue-- 
the dollar cost of forecast accuracy. That is, for a given 
demand pattern how desirable is it from an inventory cost 
perspective to pursue a more accurate forecasting model? 
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VI. COST IMPACT OF FORECAST ERROR 



The final phase of the study uses a computer simulation 
to examine how forecast error affects ordering, holding and 
stockout costs. This examination does not suggest how to 
improve forecasts but it does indicate if it is worth doing 
anything at all. 

The <Q,r> model presented in Chapter IV was a repre- 
sentation of a theoretical inventory cycle for which sta- 
tionary annual demand rates were applicable. Here the 
interest is in measuring the actual inventory costs which 
an ICP incurs when annual demand, due to forecasting 
inaccuracy, is no longer stationary and when both the 

reorder quantities (Q^) and reorder points (r^) vary for 
th 

each i cycle. To conduct such an analysis requires a 
type of "bookkeeping" computer routine capable of both 
maintaining the three inventory cost "accounts" and com- 
puting and r^ for each cycle. 

th 

For the purposes of this chapter, the 1 cycle length 
(T^) is defined as the interval of time between receipt of 
successive orders. Registering cycle order costs reduces 
to the simple task of posting the cost per order (A) once 
each cycle. The second account, holding costs, is posted 
in one of two ways. First, if no backorders occur then the 
cycle holding costs are ^T^ [ (Opening Inventory)^ + (Closing 
Inventory ) ^_^] . If a backorder does occur, however, the 
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cycle holding costs are just ^ [T.^ (Opening Inventory)^] where 
T! is that increment of for which the inventory balance 
was positive. It should be noted that in exercising this 
actual inventory cycle the opening cycle inventory may 
differ from the Q+(r-y) quantity found in the theoretical inven- 
tory cycle of Chapter IV. A more accurate representation is 
nQ^+(r^-y^) where n is an integer. This is necessary to ensure 
that the reorder point for each successive cycle will be 
reached. That is, since and r^ change with each cycle 

under certain situations an integer multiple of must be 

th 

ordered during the i cycle to enable the i+1 opening 
cycle inventory to exceed its reorder point, r^ +1 - 

As to the third account, stockout costs, two possibili- 
ties must be considered. The first possibility is that the 
backorder cost is tt per actual backorder with no penalty 
assigned for the length of time the backorder exists. 

A 

Second, the actual backorder cost, r, could also be con- 
sidered as a function of the time the backorder exists 
(T^ - T|) . Separate accounts are used to examine the effects 
of forecast error on each possibility. 

Before presenting a summary of the structure of the 
computer simulation, a measure of forecast error is required 
to demonstrate the cost effects for various magnitudes of 
forecast error. The difficulties of dealing with a single 
criterion of forecast accuracy was addressed earlier. In 
spite of the fact that it can at times be an unreliable 
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figure of merit (see Tables 3 and 4) , MSE was selected as 
the measure of accuracy. Specifically, an array of eleven 
equally spaced values ranging from 0 to 5000 was chosen 
(i.e., MSE(k) e [0,5000], k = 1,11). This range of values 
partially reflects the MSE quantities found in Tables 3 
and 4 . 

The simulation computer program (see MSESIM FORTRAN) 
borrowed the input parameters found in the INVENTOR source 
code discussed in Chapter IV (with a = 10 here) . In the 
MSESIM routine a demand value ( A^ ~ NORMAL (180 , 225) ) and a 
leadtime (LT^ ~ NORMAL ( 6 (mos) , 1) ) were generated for each 
cycle using an available random number generator (see [Ref. 
12]). Using the synthetic demand rate, the available cycle 
stock level is depleted until the cycle reorder point (r^) 
is reached. This triggers the placement of an order of size 
nQ^ such that this replenishment quantity will exceed r^ + ^. 
Stocks continue to deplete at a rate of A^ over the syn- 
thetic leadtime period (LT^) . Upon reaching the end of LT^ 
the cycle T^ is completed and the reorder quantity nQ^ 
arrives. At this point the accounts for ordering, holding 
and stockout costs are updated to reflect the costs actually 
incurred during the cycle. Next, Q^ + ^, and r^ +1 are com- 
puted using the basic algorithm found in Chapter IV except 
that A^ is deliberately falsified to reflect the forecast 
error introduced by the chosen forecasting model. (Fore- 
cast errors (FE^) were randomly generated as NORMAL (MSE (k) , 
1) deviates.) In this way, the FE^ provided by the 
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forecasting model causes a non-optimal and r i+i to 

be computed. This process continues until a time counter 
reaches a specified upper limit at which time the evolution 
is rerun for another pass. After fifty passes the next 
candidate forecast error (MSE(K+1)) is introduced and the 
entire simulation process repeats for another set of fifty 
passes. After all eleven values of MSE(K) are used the 
computer program plots the grand average for each cost 
account as a function of its corresponding value of MSE(K) 
producing the graphs which follow. 

The figures presented below are the results of the MSESIM 
computer routine just described. Again each ordinate axis 
is normalized to facilitate comparison. Also for complete- 
ness both forms of stockout costs (unit and time weighted) 
and their corresponding total variable cost graphs follow 
the ordering and holding cost plots. 

A review of Figures 5-10 shows that both ordering and 
holding costs exhibit relatively minor changes over the 
relevant range (9% decrease; 6% increase, respectively) . 

The roughly linear increases in both stockout costs are 
quite dramatic especially for time-weighted backorders where 
greater than a five fold increase is observed. Figures 8 
and 10, the two total variable cost graphs, must be viewed 
with caution for, once again, their behavior is directly 
influenced by the particular values selected for A, IC, tt , 
a . 
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NORMALIZED AVERAGE ORDERING COST 




Figure 5. 



Ordering Cost Sensitivity to MSE 
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Figure 6. Holding Cost Sensitivity to MSE 



50 



NORMALIZED AVERAGE STOCKOUT COST-UNIT 



\ 




Figure 7. Stockout Cost (Unit Weighted) Sensitivity 
to MSE 
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Figure 8. Total Cost (Unit Weighted) Sensitivity to MSE 
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Figure 9 . Stockout Cost (Time Weighted) Sensitivity 
to MSE 
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AVERAGE TOTAL VARIABLE COST-TIME WEIGHTED 




Figure 10. Total Cost (Time Weighted) Sensitivity 
to MSE 
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The conclusions gleaned from these graphical displays 
are two. First, the consequences of even severe forecast 
error on both ordering and holding costs appear relatively 
insignificant. The second and more profound finding is 
that a particular forecasting model's accuracy (as measured 
by MSE) appears to play a vital role in keeping total 
stockout costs reasonable. 
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VII. CONCLUSIONS 



A. SUMMARY AND FINDINGS 

The study began with a comparison of the SPCC forecasting 
process to an unrefined form of a newer computer-oriented 
technique. The use of several effectiveness measures showed 
neither exponential smoothing nor focus forecasting per- 
formed significantly better in forecasting actual SPCC 
demands. One outgrowth of the comparison was the prepon- 
derance of erratic-demand patterns in the ICP inventory 
population. Next a sensitivity analysis was designed to 
measure the cost effects of demand variability. This analy- 
sis indicated that reductions in demand variability gener- 
ate generous cost savings, particularly in the form of 
reduced stockout costs. Not content with viewing the erratic 
demand situation as inflexible, two research directions were 
pursued: (1) a change to the demand pattern; (2) the cost 

effects of changing the forecasting model. First the in- 
clusion of an artificial demand reporting channel produces 
considerable reductions in demand variance. A second and, 
presumably, more complicated improvement to the demand 
forecasting process would be the development of a model 
capable of accurately forecasting erratic items. To gain 
insight into the expected economic worth of any such model, 
the study concluded with an examination of the cost effects 
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of forecast error. There it was found that stockout costs 
decrease drastically as forecast accuracy improves. 

B. RECOMMENDATIONS FOR FURTHER STUDY 

Follow-on research in the areas listed below may aug- 
ment the research conducted in this study. 

1) Conduct further investigation into modifying the TIR 
reporting system by routing all consumer demands 
directly to the ICP. 

2) Investigate other causes of demand variability and 
examine the corrective procedures necessary to 
smooth input demand data. (Several possible sources 
were hypothesized in Chapter V.) 

3) Continue to analyze computer-oriented forecasting 
techniques. Since the development of the exponential 
smoot'hing concept, technological advancements in the 
computer field have lessened computational and data 
storage costs. Newer, innovative forecasting methods 
capitalize on these advancements. Research into 
dependable forecasting models which are easily under- 
stood, easily maintained and easily adapted by a multi- 
item inventory activity should receive continued 
sponsorship. The focus forecasting concept is one 
possible candidate. 
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APPENDIX A 



COMPUTER PROGRAM SOURCE CODE 
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